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Abstract 

The inclusive production rates and differential cross-sections of photons and mesons with a 
final state containing photons have been measured with the OPAL detector at LEP. The light 
mesons covered by the measurements are the n°, 77, p(770) ± , a; (782), 7/(958) and ao(980) ± . The 
particle multiplicities per hadronic Z° decay, extrapolated to the full energy range, are: 

(ra 7 ) = 20.97 ± 0.02 ± 1.15 , 

(tv>) = 9.55 ± 0.06 ± 0.75 , 

(n, n ) = 0.97 ± 0.03 ± 0.11 , 

(n p ±) = 2 AO ± 0.06 ± 0.43 , 

(nj) = 1.04 ± 0.04 ± 0.14 , 

(ry) = 0.14 ± 0.01 ± 0.02 , 

(n+) = 0.27 ± 0.04 ± 0.10 , 

where the first errors are statistical and the second systematic. In general, the results are in 
agreement with the predictions of the JETSET and HERWIG Monte Carlo models. 
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1 Introduction 



In high-energy collisions, the transition from quarks and gluons to stable hadrons can only be 
described by phenomenological models M. Among the basic features that these hadronisation 
models attempt to reproduce are the multiplicity and energy spectrum of each hadron species. 
The large sample of hadronic Z° decays collected at LEP is ideal to test these models and 
to improve their accuracy, as the initial state in this process is theoretically well understood. 
With the versatility of the LEP detectors, these measurements can be extended to most of the 



low-lying particle states [0-12 



This paper describes the OPAL measurements of the differential production cross-sections 
in hadronic Z° decays of photons and of light mesons decaying to final states containing at least 
one photon. The mesons studied are the n°, r/, p(770) ± , a; (782), 7/(958) and ao(980) ± . The 
measurements require a good understanding of both the detector response to photons and of 
the environment in which these particles are produced. For this reason, the results obtained 
with photons detected as energy deposits in the electromagnetic calorimeter and as pairs of 
tracks from photon conversions to e + e~ in the central drift chamber are first compared and 
then combined. This comprehensive study of the production of mesons decaying to photons 
leads to a better understanding of the systematic effects related to photon detection, making 
possible a reliable measurement of the inclusive production of photons in Z° decays in a wide 
energy range[|. The production cross-section of each particle is presented as a function of its 
scaled energy x E = -Eparticie/ E heam and of \og(l/x p ), where x p = Pparticie/Pbeam is the scaled 
momentum. 

Current measurements of photon and tt°, t], u and rj meson production at LEP are limited by 
experimental systematic errors [H,^-^]. Compared to these studies, the present measurements 
cover a larger fraction of the total rate of these particles, and a number of the sources of 
systematic error are different. This is the first measurement of the inclusive production of p ± 
and ag in high-energy e + e~ collisions. 

The outline of the paper is the following. The OPAL detector is briefly presented, followed 
by the description of the event selection and simulation. The following three sections describe 
the three steps in the particle reconstruction. First, photons are detected either as localised 
energy deposits in the electromagnetic calorimeter or as two tracks from a 7 — > e + e~ conversion 
within the volume of the central drift chambers. In the second step, photons are combined in 
pairs to form 7r° and 77 meson candidates. In a final step, the tt° and 77 meson candidates are 
combined with one charged track or two oppositely charged tracks to reconstruct rj — ► 7r°7r + 7r~, 
p ± — ► -R°n ± , uj — ► 7T°7r + 7r~, rj — ► r]n + n~ and a^ — ► r/^ decays. Each step is described, together 
with the corresponding rate measurements and evaluations of systematic errors. The following 
section describes the combination of the results for those particles where more than one decay 
mode is used. The resulting differential cross-sections for photons and light mesons are then 
compared to the predictions of different models and to other measurements at LEP. This is 
followed by the conclusion. 

For an inclusive measurement of prompt photons, excluding hadron decays and initial state radiation, see 
ref. El. 
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2 The OPAL detector 



The OPAL detector and its performance are described in detail elsewhere [[14]]. Only detector 
elements of importance for this analysis are described here. The central tracking system consists 
of three drift chambers which surround a silicon microvertex detector ||15|| , all within an axial 



magnetic field of 0.435 T. The silicon microvertex detector has two layers, at radii of 6.1 and 
7.5 cm from the beam axis, with an intrinsic resolution of 5 /im in the r — planeQ. A precision 
vertex drift chamber with 24 cm outer radius provides space points with a resolution of about 
50 /im in r — cf) and 1 mm in z. The jet chamber has an outer radius 185 cm and provides 
up to 159 measurements of space points per track, with a resolution in the r — <fi plane of 
about 130 /im. The resolution of the r — <fi component of the track momentum (p t ) is u Pt /pt ~ 



(0.02) 2 + (0.0015p t ) 2 , where p t is in GeV/c. In addition, charged particles can be identified by 
their specific ionisation energy loss (dE/dx) |l||. The jet chamber is surrounded by a system of 
z-chambers, thin drift chambers with a resolution of about 300 /im in the z coordinate, which 
serves to improve the determination of 9. 

The tracking detectors and the magnet coil are surrounded by electromagnetic and hadronic 
calorimeters and muon chambers. In this work, the identification of photons is performed within 
the acceptance of the barrel electromagnetic calorimeter. This consists of a cylindrical array 
of 9440 lead glass blocks of 24.6 radiation lengths thickness that covers the polar angle range 
| cos 6* | < 0.82. Each block subtends approximately 40x40 mrad 2 . The energy resolution 
is improved by correcting for losses due to showers initiated in the material in front of the 
calorimeter. These showers are detected by thin presampler gas detectors covering the front 
surface of the electromagnetic calorimeter. Time-of-flight scintillators, situated between the 
magnet coil and the presampler in the polar angle range | cos#| < 0.72 are also used to detect 
these showers. 



3 Data selection and event simulation 



This study is based on a sample of 4.1 million hadronic Z° decays collected by the OPAL 
detector at LEP at centre-of-mass energies within ±2 GeV of the Z° peak. The hadronic event 
selection [T% has an efficiency of 98.4 ± 0.4% with a background of less than 0.2%. 



The detection efficiencies for the particles under study are evaluated using 6.4 million 
hadronic Z° decays simulated using the Monte Carlo programs JETSET 7.3 and 7.4 |T^] tuned 
to reproduce the global features of hadronic events as observed at LEP |]19|,|20|. Samples gener- 



ated by the HERWIG 5.9 program [19,21] are also used for comparison. The generated events 



are passed through a full simulation of the OPAL detector |22| and are subjected to the same 
event reconstruction and selection as the data. 



2 The OPAL coordinate system is denned so that z is the coordinate parallel to the e~ beam, r is the 
coordinate normal to this axis, 9 is the polar angle with respect to z and cj) is the azimuthal angle about the 
z-axis. 
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Photons 



4.1 Photons reconstructed using the electromagnetic calorimeter 

The electromagnetic calorimeter provides the largest part of the photon sample. To resolve a 
maximum number of photons in the dense environment of hadronic jets, the location and the 
energy of the electromagnetic showers are obtained from a fit to the energy deposited in the 
individual lead glass blocks. The fit uses a description of the lateral shower profile as the sum 



of two exponentials (see for example ref. [23]) and allows a proper treatment of overlapping 
showers. The fit can also handle the cases where most of the photon energy is in a single block, 
a common occurrence for photons in the energy range from 0.1 to 0.3 GeV. 

Not all energy deposits in the electromagnetic calorimeter are due to photons. Many are due 
to ionisation or to electromagnetic and hadronic showers caused by charged particles. For this 
reason, a block lying close to the extrapolated path of a charged track is given a small weight in 
the fit, provided its energy does not significantly exceed the expectation for a hadronic shower. 

Photons may lose energy (typically about 0.2 GeV) by initiating an electromagnetic shower 
before reaching the calorimeter. The photon energy is therefore corrected using the signal 
recorded in the presampler. The efficiency for detecting these showers is further increased by 
also using the presence of signals in the time-of-flight detectors. 

A shower is retained as a photon candidate if it has at least 0.1 GeV energy in the lead glass 
calorimeter and if the energy corrected using the presampler and the time-of-flight detectors is 
at least 0.15 GeV. An acceptance cut of | cos^l < 0.75 is imposed to improve the homogeneity 
in the amount of material traversed by the photons before reaching the calorimeter. The 
momentum direction of the photon is evaluated assuming that it originates from the primary 



event vertex determined as described in 24 . 



According to the Monte Carlo simulations, the overall efficiency for photons above 0.15 GeV 
is 69% within the acceptance | cos#| < 0.75. The purity of the sample is 52%, with the most 
important background being due to energy deposits from charged particles, neutrons and K° 
mesons. In simulations of hadronic Z° decays, the photon angular resolution is approximately 
10 mrad and the energy resolution varies from 30% at 0.15 GeV to 8% at 20 GeV. 

For the rate determination, it is useful to compare results obtained from samples with 
different purities. As shown in fig. [I], the purity depends on the photon energy, the energy 
deposited in its vicinity, its distance to the closest charged track, the shower shape and the 
quality of the shower fit. The seven variables shown in fig. [I] are combined in one variable S: 

1 N 

s = AfE( 1 + ex p(( u *- c <)/*<)) , (i) 

where the index i runs over the N=7 variables v , and the parameters q and ti are chosen such 
that the power of the variable 5* to discriminate between signal and background is maximal. 
This is achieved by minimising the ratio 
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where ji s and a s are the average and rms values of S for the signal and fij, and aj, are the 



corresponding values for the background. The minimisation is performed using MINUIT [25] 
and a sample of simulated events. The contributions to S of the input variables are shown as 
inserts in figs. |l]a to 

The variable S can be interpreted as the output of a simplified artificial neural network, 
where the number of parameters optimised using a reference sample of simulated events is 
reduced to two per input variable. In this way, S is forced to depend on the global properties 
of the input variables, which should reduce problems due to the imperfection of the detector 
simulation. In the present application, the discrimination losses relative to more complex 
artificial networks are negligible. In contrast to a likelihood method, correlations between the 
input variables are partially taken into account by the simultaneous optimisation of all the q 
and ti parameters. 

The purity of the Monte Carlo photon sample, P 7 , has a smooth dependence on S which 
is easily parameterized by an analytical function. The distribution of P 1 is shown in fig. [| It 
is well reproduced by the simulation, in particular in the region P 1 > 0.5, where the signal is 
expected to dominate. The variable P 7 is used for the systematic studies of photon samples of 



varying purity described in section 4.3 



4.2 Reconstruction of photon conversions 

According to the Monte Carlo simulation, 7% of the photons in the angular range | cos 9\ < 0.75 
convert to e + e~ pairs in the volume of the central tracking chambers. It is useful to compare the 
results obtained using this sample with those obtained from photons detected in the calorimeter, 
since they are affected by different systematic uncertainties. 

The selection of photon conversions is optimised to have a high efficiency at low momentum 
and a good angular resolution, and to be insensitive to details of the detector simulation. 
Conversions are observed as two oppositely charged tracks in the central detector. The two 
tracks must have impact parameters relative to the primary event vertex greater than 300 /xm 
in the r — plane, and at least one track must have space point measurements inside the 
^-chambers. The average of the two points where the tracks are parallel in the r — <fr plane 
is taken as the point of conversion. The pair topology is required to be consistent with the 
expectation for a conversion: 

• the distance between the two points where the tracks are parallel in the r — <fi plane is 
required to be less than 1 cm, 

• the radial coordinate of the point of conversion, r, must be greater than 3 cm, 

• the reconstructed photon must have an impact parameter relative to the primary event 
vertex in the r — plane smaller than 5 cm, 

• the absolute value of the difference between r and the radial coordinate of the first space 
point measurement on either track must be less than 20 cm. 
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According to the Monte Carlo, each of these topological cuts removes less than 2% of the 
conversion sample. This loose selection is sufficient to obtain a 90% pure sample for r > 
50 cm. For r < 50 cm, the background increases because of the large track density and it 
is further required that the dE/dx measurements of the two tracks each have a probability 
greater than 1% for the electron hypothesis. According to the simulation, this cut removes 4% 
of the conversions, achieving a purity of 85% for the entire conversion sample. No further cuts 
are applied since this purity is sufficient to obtain invariant mass spectra of pairs of photon 
candidates where the background is dominated by random combinations of genuine photons. 

The direction of the photon in 9 is computed from the track parameters with the added 
constraint that the photon comes from the primary event vertex. The polar angle of the photon 
direction is required to be in the same fiducial region used for the calorimetric sample (| cos#| 
< 0.75). 

The distribution of the radial coordinate r of the conversion points is shown in figure 0. 
The Monte Carlo simulation reproduces the overall shape due to the local concentrations of the 
material in the detector. However, the numbers are not reproduced exactly: for example, in 
the data 58.0% of the conversions lie below r = 50 cm, while in the simulation this fraction is 
59.2%. Such discrepancies are considered in the following estimation of the systematic errors. 

According to the simulation, in the energy range from 0.15 to 20 GeV, the angular resolution 
on the direction of the photon conversions decreases from 10 to 1 mrad in </> and from 24 to 15 
mrad in 9, and the energy resolution is approximately constant at 4%. Approximately two thirds 
of the photon conversions in the fiducial region are reconstructed and selected, corresponding 
to an average photon efficiency of about 3%. 



4.3 Evaluation of the photon yield 

As the sizes of the photon candidate samples are large, the precision of the measured yields is 
expected to be limited by systematic uncertainties. It is therefore important to compare the 
yields derived from samples obtained with different selection procedures, using different Monte 
Carlo predictions for the photon efficiencies and different methods to subtract the background. 
Here, the size of the final systematic errors is reduced by incorporating these tests in the 
determination of the yield itself. 

In a first step, the number n(E 7 , AE y ) of photons per hadronic Z° decay in an energy bin 
of width AE^ centred at E 1 is derived using five different data samples (noted by the index i), 
three Monte Carlo samples (index j) and two methods to estimate the background (index k): 

n^( El , AE,) = N ^^ : AE ^ ~ N ^ E ^ AE y ) 

Here A^ andidates is the total number of photon candidates in the data, N^l^ is the predicted 
number of fake photons, e lJ is the efficiency for photons in that energy bin, and N z o is the 
number of Z° decays. In the Monte Carlo, the background is defined as those candidates that 
are not unambiguously associated to a single photon, or photons coming from bremsstrahlung 
radiation or decays of particles produced in interactions with the material of the detector. For 
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example, the Monte Carlo predicts that about half of the candidates above 15 GeV result from 
the overlap of the two photons from high-energy ir° decays. These unresolved photons cannot 
be counted appropriately and are therefore considered as background. 

The five data samples are the conversion sample, three calorimetric samples with require- 
ments P 1 > 0.0, 0.5 and 0.75, and the calorimetric sample with each entry weighted by P 1 . 
The three Monte Carlo samples used for evaluating and e lJ are based on the JETSET 

tunes of ref. |19| and [BO], and the HERWIG tune of ref. |19| |. The two prescriptions for the 



determination of N^l^ are to take the Monte Carlo prediction and to scale it either accord- 
ing to the number of events or the number of photons. In total, 30 energy- dependent yields 



n ,J ' fc (£' 7 , AE 7 ) are obtained. 



In a second step, the central value for n(E 7 , AE^) is obtained from a weighted average of 
the 30 analyses. The average is first performed on the background assumptions k, using the 
quadratic sums of the statistical errors on the data and on the Monte Carlo as weights. This 
results in averaged yields n JJ and the rms deviation around this mean, cr^ g , is taken as the 
systematic error on the choice of background assumptions. This error is added in quadrature 
to the total error, and the n lJ 's are averaged over the choice of Monte Carlo j, using the 
new total error as a weight. This in turn yields new averages, n l , and the error associated 
with the choice of Monte Carlo cr\ ICl which are again added in quadrature to the total error. 
Next, the same procedure is applied to the four calorimetric samples, resulting in an average 
for these samples, n cal , and a systematic error associated with the use of P 7 . In the end, two 
independent measurements of n(E 1 ,AE 1 ) are obtained, one from the calorimetric sample, and 
one from the conversion sample. At the same time, the three systematic errors associated to 
the variations of i, j and k have been calculated. The weights used to evaluate the average 
yields are also used to calculate the average error associated to each source, assuming that the 
samples being combined are fully correlated. For simplicity, this conservative assumption is 
also applied to the statistical errors on the data and Monte Carlo samples. This is justified 
because the statistical errors represent a small fraction of the total errors and, in general, there 
is a large overlap between the samples being combined. A notable exception is the combination 
of the calorimetric and conversion samples, which is discussed in section |7|. 

The averaging procedure is performed separately for each energy bin in order to take into 
account the variation with energy of the nature and the size of the systematic uncertainties. 
As a check, the order of the averaging procedure is reversed, and the resulting yields and errors 
are compared. No significant differences are observed, indicating that the different systematic 
tests are largely uncorrelated. However, the variations arising from a given test are assumed to 
be fully correlated from one energy bin to another. 

Two corrections are applied to the average yields. The first one is for the difference in 



energy calibration between the data and the Monte Carlo explained in detail in section |5.2.3 . 
The second accounts for the difference between the data and the Monte Carlo in the number of 
photons initiating a shower before reaching the calorimeter. A study of the fraction of calorime- 
ter photons with an associated signal in either the presampler or time-of-flight reveals that the 
Monte Carlo underestimates by 2% the probability of initiating a shower before reaching the 
calorimeter. The efficiency is corrected for this effect. It is important only for photons with an 
energy comparable to the average energy lost before reaching the calorimeter, i.e., about 200 
MeV. 
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The numbers of photons per hadronic Z° decay in the energy range xe > 0.003 obtained 
from the calorimetric (7) and conversion (7J measurements are given in table [3] together with 
the values of each systematic uncertainty. The sources of these uncertainties are: 



• The statistical error on the Monte Carlo samples used to calculate the efficiency 

• The variations observed using different Monte Carlo samples, obtained from the averaging 
procedure. 

• The error associated with P 7 , obtained from the averaging procedure. 

• The variations observed when using the different background assumptions, obtained from 
the averaging procedure. 

• Some of the background comes from photons produced in interactions with the material 
of the detector. There are indications that this effect is not exactly reproduced by the 
Monte Carlo; see for example ref. |2| and the discussion on the electron bremsstrahlung 
in section |5.2.3| . Therefore, half of the Monte Carlo prediction for this background is 
conservatively taken as an uncertainty. 

• As the measured yields depend on the exact energy calibration, the analyses are repeated 
by shifting the energy scale by 1% and the difference in the rates is taken as a systematic 
error. The size of this shift is justified by the calibration studies described in the following 
section on tt° and 77 reconstruction. 



The variation of the cuts on P 7 does not cover all sources of systematic uncertainties on the 
quality of the simulation. The following systematic errors are considered: 

• For the calorimeter data, the error labelled simulation in table [l] is the quadratic sum of 
the uncertainty on the correction of the probability of initiating a shower before reaching 
the calorimeter and the yield variation observed when the criteria for the association of 
charged tracks to energy deposits are varied. 

• For the conversion data, the error labelled simulation is obtained by removing each of the 
selection cuts in turn. The differences between the number of accepted tracks in the data 
and the Monte Carlo obtained when each cut is removed are added in quadrature. 



For the calorimetric measurement, the largest source of systematic uncertainty is the dif- 
ference between the efficiencies derived from the three simulations (table [I]). In contrast to 
ref. I [L 0|| , where the most important error arises from the choice of using either JETSET or 
HERWIG for the determination of the efficiencies, here the largest difference is observed be- 
tween the two samples generated using the JETSET versions tuned in refs. [^9| and pO . The 



efficiencies obtained with the HERWIG 5.9 sample are consistent those obtained with the first 
JETSET sample. Notable differences between the two JETSET samples are the inclusion of 
L = 1 mesons |n| and consequent changes to the uo and 7/ rates by factors of 1.5 and 4.5, re- 
spectively. These differences affect the photon efficiencies due to the presence of neighbouring 
particles. Another difference between the two samples is the version of the detector simulation 
program [p2|. 
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For the conversions, the largest source of uncertainty is that associated with the simulation 
of the selection cuts (table [I]). The error is 6.8% over most of the energy range and is dominated 
by effects related to the inadequacy in describing the distribution of the radial coordinate of 
the conversions. At the lowest energy it increases to 10% due to uncertainties in reconstructing 
tracks with small curvature radii. 

As a consistency check, the photon yields obtained from the calorimetric and the conversion 
sample have been compared in 10 data sub-samples corresponding to different data-taking 
periods, spanning six years of operation of OPAL at Z° energies. The calorimeter efficiency 
varies by less than 1% and the number of conversions varies with an rms of 2.6%, well within 
the systematic errors estimated for those channels. 

The differential cross sections of inclusive photon production as a function of xe and 
log(l/x p ) are presented, interpreted and discussed together with those for the light mesons 
in section |S|. 



5 The channels 7r° — ► 77 and rj — ► 77 

In this section, 7r° and r\ mesons are reconstructed as pairs of photons. The branching ratios of 
the decays vr° -> 77 and 77 -> 77 are (98.80 ± 0.03)% and (39.25 ± 0.31)%, respectively f2§. 
The numbers of it and rj mesons in the data and Monte Carlo samples are determined from 
fits to the invariant mass spectra of the photon pairs. The selection of candidates is presented, 
followed by a description of the fits to the invariant mass spectra and the determination of the 
meson yields. As in the photon case, the yields and some systematic errors are obtained by 
averaging results based on different samples and using various analysis methods. 



5.1 7T° — > 77 and 77 — > 77 selection 

The 7T° and rj candidates are obtained by combining in turn all pairs of photon candidates. In 
the rj selection, the energy of each photon is required to be larger than 0.3 GeV. At this stage, 
the combinatorial background is large, and it is not possible to extract the ir° and r\ yields from 
the invariant mass distribution of the photon pairs. 

To reduce the background, the probability P 77 that a photon pair comes from a 7r° — > 77 
decay as a function of a set of input variables is estimated using the same method as for P 7 
in section |4.1| . The input variables for each photon are those shown in fig. |l|a, b, d, f, and g, 



together with the opening angle 0^ of the photon pair, the number of additional photons in 
cones of opening angle around each of the two photons and the helicity angle of the photon 
calculated in the 77 rest frame. The invariant mass of the pair is excluded from the input 
variables so that it can be used later on for extracting the yields. The same method is applied 
to the decays involving conversions^ and the three functions PZl ~* 77 , PL ~* 77c and ~* 7c7c 
are determined separately. Similarly, the equivalent functions for the rj, P^ 7- * 77 , P 7 7 ~~ >77c and 



3 For the conversions, the input variables corresponding to fig. ^b, d, f, and g are not relevant and are therefore 
not used. 
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P 77 ~^ 7c7c , are evaluated separately. The purity of the 7r° Monte Carlo sample is evaluated using 
a Am— 0.1 GeV/c 2 -wide invariant mass window centred at mo =0.135 GeV/c 2 . For the rj, the 
purity of the Monte Carlo sample is evaluated using Am=0.2 GeV/c 2 and mo=0.5475 GeV/c 2 . 

The probability that a pair of photons i,j is part of the 7r° or 77 signal depends not only on 
the P 77 value for this pair, but also on the signal probability for any other pair involving either 
i or j. Assuming that P 77 is indeed equal to the signal probability, and assuming that the P 77 
values of all photon pairs in an event are uncorrelated, the probability that two photons % and 
j are not related to any other photon in the event is: 

Kr, = II [1 - G ik (m ik )} x [1 " G 3l {m 3l )} . (4) 
Here is the invariant mass of the pair of photons i and k, and 

G ik (m ik ) = ^3 (5) 

is the P 77 variable modified to take into account the mass dependence of the purity of the ir° 
and T) signals, assuming that the invariant mass peaks are normal distributions of width a, 
which are determined from the Monte Carlo. The products in eq. |] are performed over all pairs 
retained by either the 7r° or the 77 selection. If a pair passes both selections, the larger of the 
7T° or rj signal probability is retained. The combined probability that the photons i and j are 
from the same 7r° or 77 decay and that they do not take part in any other decay is taken as 
P 7i7j = P 7i7j x P 7 . 7j . According to the simulations, cutting on P 77 instead of P 77 reduces the 
combinatorial background by approximately 10% in the case of the 7r° and by as much as 30% 
in the case of the 77. 

Since the 7r° and rj yields are determined from fits to the invariant mass spectra, the effect 
of a cut on P 77 on the shape of the invariant mass distributions has to be studied. The value of 
P 7i7j has a monotonic correlation with the invariant mass of the pair m^, which arises primarily 
from its dependence on the opening angle % and the number of additional photons in the cones 
of opening angle Oij around i and j. The shifts in the position of the 7r° and rj mass peaks 
induced by cuts on P 77 are of the order of a few MeV/c 2 and are well reproduced by the Monte 
Carlo. More importantly, it has been verified in the Monte Carlo that the cuts on P 77 do not 
produce a fake 7r° or 77 peak in the invariant mass distribution of the combinatorial background. 

Fig. I shows the distributions of P 77 for all channels and indicates that their shapes are well 
reproduced by the Monte Carlo. A more detailed discussion of possible differences will follow 
in section |5.3| . The Monte Carlo predictions for the efficiency for detecting 7r° and 77 mesons 
using cuts of 0.1 and 0.05 on P 77 , respectively, are shown in fig. |5|. The efficiencies include all 
effects, including the known branching ratios p6| , so that the yields obtained from the 77, 77 c 
and 7 C 7 C samples can be directly compared. 



5.2 Analysis of the invariant mass spectra of photon pairs 

As in the photon analysis, the determination of the 7r° and rj yields is repeated using different 
data samples and analysis procedures, and the results are averaged to obtain the central values. 
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Therefore, the fits to the invariant mass spectra are performed separately for the ir° — > 77, 
7T° — > 77c, tt° — > Jdc, r] — > 77 and 77 — > 77 c samples. The channel r\ — > 7 C 7 C is not used because 
of low statistics. The fit is systematically repeated using five values of the cut on P 77 , two 
parameterizations of the background and two parameterizations of the signal peaks. In addition, 
three Monte Carlo samples are used for the determination of the n° and rj efficiencies and, for 
one of the background parameterization methods, the background shape. These variations are 
described below. They amount to a total of 60 different methods to extract the rates for each 
channel. 

5.2.1 Cut variations 

The cut on P 77 is varied among the values 0.1, 0.2, 0.3, 0.4 and 0.5 for the 7r° analysis, and 
0.05, 0.10, 0.20, 0.30 and 0.40 for the rj analysis. The variations of the P 77 cut are chosen so 
as to result in a change in acceptance of at least a factor two, in order to test how well G 77 (m) 
is simulated by the Monte Carlo, and also to provide a wide variety of background shapes and 
levels. 

5.2.2 Background parameterization 

Figs. H and |7] show the invariant mass spectra for all five channels. The two background 
estimations are also shown. In the first one, the background shape is taken from a simulation 
and is normalised to the number of counts outside the signal region. Possible differences between 
the shape in the data and in the Monte Carlo are taken into account by adding a linear 
background to the fit. In the second method, the background is fitted using a second-order 
polynomial. An additional source of background is considered in the region just above the 77 
peak. The reflection from u — > 771" decays are taken into account using a Gaussian with a 
normalisation allowed to vary and a mass and width of 730 and 80 MeV/c 2 , respectively. This 
peak is caused by the kinematic correlation between one of the photons from the 7r° and the 
direct photon from the u, and cannot be absorbed by the polynomial background. 

The two background parameterizations are complementary; the first one takes into account 
all the features of the background shape predicted by the simulations, while the second does not 
depend on the details of the Monte Carlo. The area of the signal peaks obtained with the two 
methods are in general not identical, as shown in figs. |6| and [7|. However both methods should 
give the same yields if the efficiency is determined by applying the same parameterization to 
the Monte Carlo. As a result, the errors due to the background parameterization are smaller 
than is suggested by the difference between the two background estimates in figs. || and [7|. 
Part of this difference comes from an ambiguity in the definition of the background in the 
Monte Carlo. In a few percent of the shower reconstructed in the calorimeter has 

contributions from more than one incident particle and cannot be clearly associated to any one 
of them. However, the contribution of a photon to this shower might be important enough 
that when it is combined with the photon coming from the same 7r° decay, the invariant mass 
of the pair may be very close to the mass of the it . Even though such pairs may produce a 
small 7T° peak, the first method consider them as being part of the background while the second 
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will tend to treat them as part of the signal. The averaging procedure takes into account this 
ambiguity and the related uncertainty. 



5.2.3 Signal parameterization 

Fig. H shows the shape of the mass peak for the samples 7r° — ► 77, n° — ► 77 c , 7r° — > 7 C 7 C and 
77 — > 77 in the data and in the simulation. In all cases the distributions are obtained by sub- 
tracting the fitted background from the raw spectrum. In the fits, the shapes are alternatively 
parameterized either as a Gaussian function, or as the more complex functions described below. 

With high statistics samples such as those in fig. ^ the non-Gaussian structure of the signal 
peaks is apparent. The peaks for the n° — ► 77 and 77 — > 77 decays are each described by two 
Gaussians centred at the same mass. The centroid, widths and amplitudes of the two Gaussians 
are determined independently for the data and the Monte Carlo. The centroids and widths are 
well described by the simulation (fig. [|a and d). In the case of 7r° — > 77 c and 7r° — > 7 C 7 C decays 
(fig. P|b and c), bremsstrahlung from the conversion electrons produces a pronounced tail toward 



low invariant masses |p7 |. In the fit, this tail is described by an exponential convoluted with the 
fitted mass resolution. The amplitude and decay constant are determined independently for 
the data and the Monte Carlo. In the case of the 77 — > 77 c signal, the tail is neglected because 
of the low statistics. 

The comparison of the centroids of the 7r° and r\ signals in the data with those obtained 
in the simulations provides an important check of the energy calibration. The corrections for 
all five channels are of order 1% over most of the energy range. This sets the scale for the 
uncertainty on the energy calibration. 

The quality of the detector simulation and its impact on the analysis can be assessed by 
comparing the widths and the tails of the mass peaks in the data and in the Monte Carlo. The 
Monte Carlo reproduces well the peak shapes (fig. |8|), except for the tail toward low masses 
for combinations involving one or two low-energy conversions. Such tails have been noticed 



before p7| an d they are due mostly to bremsstrahlung of the conversion electrons. The Monte 
Carlo simulates correctly the slope of these tails, but underestimates their amplitude by as 
much as a factor of two. The analysis should not be affected by an inadequate modelling of 
these tails, since they are free parameters of the fit. In addition, the systematic error includes 
the effect of neglecting them by assuming Gaussian peak shapes, and the effect is further tested 
in section |7| by comparing the results derived from the 7 and 7 C samples, and from the 77, 77 c 
and 7 C 7 C samples. 



5.2.4 Monte Carlo simulations 



The three Monte Carlo samples described in section |4l| are alternatively used to evaluate the 
7T° and rj efficiencies. When the first background parameterization method is used, the shape 
of background is taken from the corresponding simulation for consistency. 
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5.3 Determination of the 7r° and 77 yields 

The 7T° and 77 yields and their errors are determined separately for the channels 7r° — ► 77, 
7T° — ► 77c, 7T° — > 7c7c, 77 — > 77 and 77 — > 77 c using the averaging method described in section O . 



As the numbers of mesons -Candidates are obtained from fits to the invariant mass spectra, the 
background iVbkg entering eq. |3| is the contribution from 7r° produced in interactions with the 



material of the detector, which is taken from the simulation. As in section [4.3| , the yields are 
corrected for the difference between the Monte Carlo and the data in the energy calibration 
and in the probability that a photon initiates a shower before reaching the calorimeter. 

The numbers of 7r° and 77 mesons per hadronic Z° decay in the energy ranges covered by 
the present measurement are given in tables ^] and [3] together with the values of each of the 
systematic uncertainties. These are: 

• The statistical error on the Monte Carlo samples used to calculate the efficiency. 

• The variations observed in the averaging procedure using different Monte Carlo samples. 

• The error associated with P 77 , obtained from the averaging procedure. 

• The variations observed in the averaging procedure when using the different background 
parameterizations. 

• The variations observed in the averaging procedure when using the different signal pa- 
rameterizations. 

• The uncertainty associated with the tt° produced in interactions with the material of the 



detector, evaluated as in section 4.3 



The uncertainty associated with the energy calibration, evaluated as in section PO 



• The first method to evaluate the background requires the definition of a mass range used 
for the normalisation of the Monte Carlo prediction to the data. This procedure is not 
exact because of the presence of tails in the invariant mass distributions of the signal. 
The bias on the yields resulting from the choice of the normalisation range is estimated 
using the Monte Carlo, and its size is taken as the systematic error associated with this 
method of evaluating the background. 

• The simulation uncertainties not covered by the variations of the cut on P 77 , i.e., the errors 
labelled simulation in table [II propagated to the 77, 77 c and 7 C 7 C samples according to 
the number of calorimetric and conversion photons in the pairs. 

As a consistency check, the position of the tt° and 77 mass peaks has been measured in 10 
data sub-samples corresponding to different data-taking periods, spanning six years of operation 
of OPAL at Z° energies. The energy scale varies by less than the 1% systematic uncertainty 
ascribed to it. The fluctuations in the extracted n° and 77 yields are of the same size as those 



observed in section [4.3| for the numbers of calorimetric and conversion photons. 



The dominant systematic uncertainties on the ir° — > 77, tt° — > 77 c and 7r° — > 7 C 7 C yields 
(table 0) are those that also affect the corresponding 7 and 7 C measurements (table [I]), namely 
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the variation of the efficiency in the different Monte Carlos for the calorimetric data, and the 
simulation uncertainty associated with reconstruction cuts for the conversion. However, the 
size of these uncertainties are approximately the same as those associated with the fits to the 
invariant mass distributions. 

The largest systematic error for the rj —>■ 77 channel comes from the variation of the cut on 
P 77 (table |3|). This is due to the difference between the shape of the P 77 distributions in the data 
and in the Monte Carlo, which becomes clearly visible above P 77 > 0.5 in fig f|b. Several checks 
on the 77 sample were performed in order to understand the origin of this discrepancy. A possible 
explanation is that the Monte Carlo simulations underestimate the number of isolated 77 mesons. 
This was verified with 77 samples selected requiring that there should be no other photons or 
charged tracks within a cone of half-angle 15° from the 77. The simulation underestimates the 
number of these isolated 77 mesons with xe > 0.1 by a factor 2.07 ±0.11, while the combinatorial 
background is well reproduced. The same problem affects tt° mesons, for which this factor is 
measured to be 1.99 ± 0.05. The factors are the same for 77 and 77,- samples. They are similar 
to those observed in ref. |p8|| . Thus, there is clear evidence that the Monte Carlo simulation 
underestimates significantly the number of very isolated 7r° and rj mesons. These mesons are 
systematically associated with large values of P 77 . Fortunately, they constitute only a few 
percent of the total sample of identified mesons and their impact on the inclusive rate is very 
small. However, they are an indication of yet unexplained shortcomings in current Monte Carlo 
models, and have to be taken into account in the evaluation of the systematic error. In the 
present analysis, the error associated with the variation of P 77 corresponds to a decrease of 
acceptance by a factor of two and probes the behaviour of mesons with significantly different 
environments. An additional check is performed in sect by comparing the 77 and 7r°7r + 7r~ 
channels which have different sensitivity to variations of the cut on P 77 . 

The differential cross sections as a function of xe and log(l/a;p) are presented, interpreted 
and discussed together with those for the photons and the other mesons in section [S[ 



6 The decay channels 7r ± 7r , 7r ± ?7, 7r + 7r 7r and 7r + 7r 77 

The reconstruction of tt° and 77 mesons offer the possibility to reconstruct the dominant — > 
tt°-k ± and a — > tjtt^ decays, the r/,uj — > -k tt + it~ decays with branching ratios of (23.2 ± 0.5)% 
and (88.8 ± 0.7)% and the 7/ — > 7]tt + tt" decay with a branching ratio of (43.7 ± 1.5)% ||26|| . 
The method to evaluate the meson yields follows closely that used in section Namely, the 
numbers of mesons are obtained from fits to the invariant mass spectra of the meson decay 
products and the final results and some systematic errors are calculated by averaging the yields 
obtained using different data samples and various analysis methods. Two notable changes with 
respect to section are the slight modification of the 7r° and 77 selection described below and 
the adaptation of the fit procedure to the properties of the observed signals and backgrounds 
presented in sections |6]^ to |6~5[ 

The decay channels tt^tt , tt 77, tt + tt~7t° and 7i + 7i^i] are reconstructed by combining the 
7T° and 77 candidates with either one charged track or two oppositely charged tracks. The 
charged tracks are required to have at least 40 measured space-points in the jet chamber, 
a momentum component perpendicular to the beam axis of at least 0.15 GeV/c, an impact 
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parameter relative to the primary event vertex of less than 0.5 cm in the r — plane and 20 
cm along the z direction. In addition, the dE/dx measurement must have a probability greater 
than 1% for the pion hypothesis. The 7r° and r] selections are improved by using 

G 77 (m) = = ^ (6) 

^7 + (l-^7)^exp(i(^)2) 

instead of -P 7i7j as the selection variable, using the same values of mo, u and Am as in eq. |^. 
The 77, 77 c and 7 C 7 C samples are simply summed because their purities are approximately the 
same for a given value of G 77 (m). The resolution on the ir° and rj momenta and energies are 



improved by constraining their invariant masses to their nominal values [26 using a kinematic 
fit. The four- momenta of the charged and neutral particles are added, and the invariant mass 
of the system is evaluated. 



6.1 Determination of the particle yields 

As in the previous section, the determination of the meson yields is repeated using different 
data samples and analysis procedures and the results are averaged to obtain the central values. 
Similarly, the fits are systematically repeated using different values of cuts on G 77 (m), different 
parameterizations of both the background and the signal and using different Monte Carlo 
samples for the determination of the efficiency and, in some cases, the shape of the background. 
For channels involving a 7r°, the cut on G 77 (m) is varied among the values 0.1, 0.2, 0.3, 0.5 
and 0.7. For channels involving an rj in the final state, the cut values are 0.05, 0.1, 0.2, 0.3 
and 0.4. Since the largest deviations in the predicted efficiencies observed in the previous two 
sections are those obtained when comparing the two JETSET samples, only these two Monte 
Carlo samples are used in this section. In the following, the variations in the parameterization 
of the signals and the backgrounds are presented channel by channel. 



6.2 Analysis of the 7r°7r invariant mass spectra 

Fig. [5] shows the invariant mass spectra of 7r°7r ± combinations for the entire energy range for two 
different values of the cut on G 77 (m). The extraction of the p yield from the Tr ^ invariant 
mass spectra is complicated by the large width of the resonance, by uncertainties regarding 
its exact shape and by the reflection from u — > 7r°7r + 7r _ decays. An additional complication 
is the presence of partially reconstructed p£ decays, where only the charged pion and one of 
the photons come from the decay of a p meson while the other photon candidate has another 
origin. It is therefore particularly useful to compare yields obtained using different cuts on 
G 77 (m) since this variation produces significant shifts in the position of the maximum of the 
background shape (fig. |9|) and also changes the relative number of partially reconstructed p ± 
decays. 

Two methods are used to evaluate the background. In the first, the background shape is 
taken from a simulation. It is normalised to the number of counts outside the signal region. 
This method is used to obtain background-subtracted invariant mass spectra. In the second 
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method, the background is parameterized as: 

f{m) = Pl {A m ) p * x exp{p 3 A m + p 4 A 2 J , (7) 

where A m = m — m n o — m^± , m is the invariant mass of the 7r°7r ± system and the parameters 
Pi to P4 are determined in the fits to the data. A Gaussian representing the reflection from 
uo — > 7r°7r + 7r~ decays is added to this shape, with a width fixed to the Monte Carlo prediction. 
The amplitude and centroid of the Gaussian are free parameters, in order to absorb possible 
imperfections in the modelling of the background in this region close to the tt ^ threshold. 
The simulation predicts that the reflections from 77 — > 7r°7r + 7r~ and K*(892) ± — > K ± ir° decays 
are small; they are not included in the fit. The background shapes obtained using the two 
methods are shown in fig. |9|. Also shown are the data before and after the subtraction of the 
average of the two background shapes. The p ± resonance is clearly observed, albeit at a slightly 
lower mass than in the Monte Carlo. 

The shape of the p resonance may be more complex than a relativistic Breit-Wigner. This 
is the case for the p° meson &E9|, where the observed deviations from this shape may be due 



to residual Bose-Einstein correlations, to interference between the amplitudes of the p and the 



7T7T background and to interference with the uo [30]. The first two effects should also affect the 
p . The most apparent sign of this distortion is a shift towards low mass of the maximum of 
the resonance. For this reason the position of the pole of the resonance is a free parameter in 



the fits. Following ref . |30[ , the modification of the p shape is taken into account by multiplying 



the relativistic Breit-Wigner by a factor 

I(m,C) = 1 + c^^, (8) 

ml 

where the parameter C is to be determined from the data. The values of C obtained in fits to 
the data in different energy bins vary between 0.2 and 1.0, corresponding to shifts of —10 to —30 
MeV/c 2 in the position of the maximum of the resonance and consistent with the observations 
of ref. [0 for the p°. 

The p ± yields are evaluated either from the results of fits to the invariant mass spectra 
or by numerical integration of background-subtracted spectra. In the fits, the resonance is 
parameterized as a relativistic Breit-Wigner convoluted with the experimental mass resolution 
and multiplied by the factor J(m, C). The fits are repeated with the parameter C being either 
fixed to zero or left as a free parameter. In addition, the experimental resolution is either fixed 
to the Monte Carlo prediction or left as a free parameter. These variations affect the number of 
mesons contained in the tails of the resonance. To address this problem in a consistent way, the 
numbers of p ± mesons are obtained by integrating the fitted resonance shapes over the same 
range used for the integration of the background-subtracted spectra, from 0.39 to 1.15 GeV/c 2 . 
All numbers in the data are then multiplied by 1.08 ± 0.03, which is the average correction 
needed to extrapolate the fitted shapes to the range from the ti ti ± threshold up to 1.5 GeV/c 2 . 
The integration range corresponds to the mass range to which the resonance is truncated in 
JETSET and therefore the yields extracted from the simulations are not multiplied by this 
factor. In addition, JETSET uses a non-relativistic Breit-Wigner shape, without any of the 
correlations or interference effects just discussed and for this reason the shape of the resonance 
is determined separately for the data and the Monte Carlo samples. 

Another important consideration in the fits to the invariant mass spectra is the treatment 
of partially reconstructed p mesons. If a photon related to the p decay contains most of the 
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energy of the ir° and is combined with any low-energy photon, the invariant mass distribution 
of the system will form a broad peak under the signal. The height of this "bump" is correlated 
with the p ± yield, but it is wider than the peak of fully reconstructed mesons. The influence 
of partially reconstructed p mesons is taken into account in the analysis and in the systematic 
errors in the following ways. The bias in the extracted yields is evaluated by comparing the 
results of fits to the invariant mass spectra of Monte Carlo samples where the partially recon- 
structed p ± are included or not. In addition, the simulation predicts that the importance of 
partial reconstruction decreases as the cut on G 77 (m) (and therefore the quality of the selected 
photons) increases. The comparison of the rates measured with different cut values is thus an 
additional test of how well this effect is simulated. 



6.3 Analysis of the 7/71^ invariant mass spectra 



Fig. [TCj shows the invariant mass spectra of rjir^ combinations for the entire energy range for 
two values of the cut on G 77 (m). The analyses of the r/n^ 1 and 7r°7r ± invariant mass spectra 
are similar. The same parameterization for the combinatorial background is used (eq. |7|) and 
a Gaussian is added to represent the reflection from 7/ — > i]tt + tt~ decays. However, because of 
the low statistics the description in the fit of the a^ resonance is not as detailed as for the p ± . 
The ag peak is parameterized by a simple Gaussian. The fitted mean and width of the signal 
peak are 990 ± 12 and 51 ± 9 MeV/c 2 , respectively, in agreement with the nominal mass of 



983.5 ± 0.9 MeV/c [p^] , and consistent with the expected range of 50 to 90 MeV/c for the 



width of the resonance in the 7771" channel |32[ . To control possible biases due to the Gaussian 
assumption, the fitted yields are used in the averaging procedure together with the integrals of 
background-subtracted spectra in the range from 880 to 1120 MeV/c 2 , where the background 
shape is either taken from the simulation or from a fit of the analytical background shape to 
the data outside the signal region. 



In fig. [TO, the data are shown together with two background distributions obtained either 



by the fit to the data or by scaling the Monte Carlo prediction. Also shown is the data after the 
average of the two background shapes is subtracted. A peak is observed, with a position and a 
width comparable to the prediction of the Monte Carlo for the a^ signal. In the JETSET Monte 
Carlo the a^ resonance is a non-relativistic Breit-Wigner with a pole at m = 983 MeV/c 2 and 
a width r = 57 MeV/c 2 , the distribution being truncated at m ± 50 MeV/c 2 . Due to the 
severe truncation and taking into account the experimental mass resolution, the simulated signal 
distribution shown in fig. [TCj resembles closely that of a Gaussian of width a ~ 32 MeV/c 2 . 



According to the simulation, the effect of partially reconstructed 77 mesons is negligible. This is 
because the 77 selection imposes more stringent requirements on the photon quality than does 
the 7T° selection. 



6.4 Analysis of the tt tt+tt invariant mass spectra 

Fig. [TJJa and c show the invariant mass spectra of 7r°7r + 7r~ combinations, for the entire energy 
range, in the region close to the 77 and u signals. Compared to the invariant mass spectra, 
the 7r°7r + 7r~ analysis benefits from the narrow widths of the 77 and uj states. However, it suffers 
from the reduced meson rates and the additional combinatorial background. 
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In the fits to the invariant mass distributions, the peaks from tt°7t + tt~ decays of the r\ and 
the uj are each well reproduced by two Gaussians sharing the same centroid, their relative widths 
and areas being determined from the simulation. In the case of the uj, this double-Gaussian 
also helps to account for the tails of the Breit-Wigner distribution of the resonance. The fitted 
parameters of the peak are the area, its centroid and its rms width. 

The combinatorial backgrounds for the rj and the uj are described by a second- and a 
third-order polynomial, respectively. The shape of the fit component representing partially 
reconstructed mesons differs from the one used for the p ± . The rj and uj momenta are shared 
amongst three daughters instead of two as in the case of the p . Consequently, the cases where 
all particles except one of the two photons come from the same decay result in a bump with 
a narrower width, lying systematically at masses above the peak of fully-reconstructed decays. 
The adopted shape is an exponential tail at high masses convoluted with a Gaussian with a 
width equal to that of the signal peak. The exponent of the tail and its area relative to the 
signal peak are taken from the Monte Carlo prediction. 

Two tests are performed to verify that the bump due to partially reconstructed decays is 
properly described by the Monte Carlo. The first is to compare the rates obtained with different 
cuts on G 77 (m). As the tt° probability increases, the size of the bump relative to the signal 
decreases. The second test uses the characteristic matrix elements of the r] and uj decays to 
produce invariant mass spectra with almost no background. Within an experimental precision 
of a few percent p4j , the decay transition probabilities and A^ are proportional tof]: 



K « 1--^- (9) 

J 0,max 

Ac ex \p!*pl\\ (10) 

where T * is the kinetic energy of the tt° in the 7r°7r + 7r~ rest frame, T * maa , is its maximum possible 
value and p+ (p*) is the momentum of the positively (negatively) charged pion in the 7r°7r + 7r~ 
rest frame. A^ and A^ are normalised such that they vary from to 1. Random combinations 
of three pions distributed according to phase space result in a flat A distribution. Therefore 
the ratio of the A distributions for the signal and the combinatorial background should be 
proportional to A. Extracting from the invariant mass spectra the component proportional 
to A, one obtains distributions with the combinatorial background subtracted, and where the 
shape of the signal and of the bump of partially reconstructed mesons can be evaluated with 
more precision. 

The method to extract the component proportional to A in the data does not depend on its 
modelling in the Monte Carlo. The A distribution of the candidates is sampled as a function of 
their invariant mass, m. The behaviour of phase space, F(X,m), including acceptance effects, 
is taken from combinations with invariant masses above and below the mass peak. In the 
signal region, it is interpolated using a polynomial function. In each invariant mass bin, the 
data are fitted with two components, one proportional to the phase space behaviour, F(A,m), 
and the other proportional to XF(X,m). Figs. |TT| b and d show the components proportional 
to XF(X,m) extracted from the data in figs. |TT| a and c, respectively. The distribution for the 
uj signal is slightly asymmetric due to a small excess at high mass which is explained in the 
Monte Carlo as the bump of partially reconstructed mesons (dashed line in fig. |TT|b). Tests 

4 As formulated here, must be multiplied by the Tq phase-space dependence. 
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performed on the simulations show that the components proportional to XF(X, m) are excellent 
approximations of the signal distributions. In both data and Monte Carlo, the meson yields 
obtained from fits to the total invariant mass distribution or to the component proportional to 
A agree within 3%. 

In JETSET, the matrix element of the decay rj — > 7r°7r + 7r~ is not simulated. To reproduce 
the data, the signal events in the Monte Carlo are weighted by a factor^ of 2X V . Due to the 
strong correlation of A^ with the 7r° energy, this correction changes the total efficiency for the 
detection of r] — > 7r°7r + 7r~ decays by as much as 20%. 



6.5 Analysis of the r]7T + 7T invariant mass spectra 



Fits to the 777r + 7r~ invariant mass spectra for the entire energy range are shown in fig. [12|, for 
the data and the Monte Carlo. Given the small statistics, a Gaussian is found to describe 
adequately the r( signal. As for the a^, the contribution from partially reconstructed r\ decays 
does not need to be parameterized in the fit. In contrast to the three-body decays of the r\ 
and u, the Dalitz plot for the decay rf — > i]7i + 7i^ is closer to phase space [f3H] and the matrix 
element of the decay cannot be used to obtain background-free distributions. 



6.6 Determination of the meson yields 

The p ± , aj, 77, u and 77' yields and their systematic errors are determined using the same 
averaging method as in sections |4.3j and |5.3| . The yields are corrected for the known branching 
ratios of the different decay modes |[26|| , except for ao — > r/^, where a branching ratio of 
90 ± 10 % is assumed 



The numbers of mesons per hadronic Z° decay in the energy ranges covered by the present 
measurement are given in table ^ together with the values of each systematic uncertainty. These 
are: 



• The statistical error on the Monte Carlo samples used to calculate the efficiency. 

• The variations observed using different Monte Carlo samples, obtained from the averaging 
procedure. 

• The error associated with G 77 (m), obtained from the averaging procedure. 

• The variations observed when using the different background parameterizations, obtained 
from the averaging procedure. In the case of 7r°7r + 7r - decays, this includes the variations 
observed when the matrix elements A^ and A^ are used. 

• The variations observed when the mass resolution is either fixed to the Monte Carlo 
prediction or left as a free parameter. In the case of the p ± , this also includes the 
variations observed when C is either fixed to zero or left as a free parameter. 

5 The normalisation is chosen such that the integral of the signal, J* F(X, m)d\, is 1. 
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• For the p , the error associated with the extrapolation of the resonance beyond the mass 
range described in section |6~^ . 

• The simulation uncertainties not covered by the variations of the cut on G 77 (m). The 
contribution of the 7r° and 77 are the simulation errors quoted in tables ||] and |^. These are 
added in quadrature and correspond to a 2% uncertainty on the charged pion selection. 

• The uncertainty due to partially reconstructed mesons, evaluated with the Monte Carlo 
by comparing the results of fits to mass distributions with and without their contribution. 
The yield difference is taken systematic error. 

• In the case of the decays 77 — > 7r°7r + 7r~, 7/ — > r]n + n~ and a^— > rjir^, the uncertainty on 
the branching ratio is greater than 1%, and is therefore included. 

The differential cross sections as a function of xe and \og{l/x p ) are presented, interpreted 
and discussed together with those for the photons and the other mesons in section [8[ 



7 Combination of channels 

The calorimeter and conversion data are compared. The ratio of the total photon rates ob- 
tained using identified photon conversions and the calorimeter data is 1.010 ± 0.002 (stat.), in 
agreement well within the size of the systematic errors on the two samples (4.7% and 8.3%; 
see table [l]). Comparison of the 7r° rates measured in the energy range common to the three 
channels 77, 77 c and 7 C 7 C also yields ratios consistent with one. The ratios of the rates 77 c /77 
and 7 c 7 c /77 are 0.97 ± 0.02 (stat.) and 0.96 ± 0.05 (stat.), respectively, in perfect agreement 
given the estimated 8.4%, 12.2% and 23.6% errors on the rates for the 77, 77 c and 7 C 7 C chan- 
nels, respectively (table |2|). The agreement is also excellent for the 77 c and 77 channels for the 
77, despite a 34.3% uncertainty on the 77c channel (table [|). 

Given the positive result of these consistency tests, a weighted average of the differential 
cross-sections measured using the different channels is performed. The weights are taken as 
the inverse of the square of the total errors. The systematic errors on the combined results 
are obtained by assuming that the individual systematic errors are completely correlated. In 
the few cases where the individual measurements and their errors are more than one standard 
deviation away from the averaged result, the error on the average is scaled using the method 
of ref. 

The 77 differential cross section measurements based on the 77(77 c ) and 7r°7r + 7r~ channels 
have comparable total errors (13.9% and 12.9%, respectively; see tables |3] and |j) and the sys- 
tematic uncertainties are largely uncorrelated. The two measurements are combined assuming 
that the errors associated to the track and cluster simulation and to the variations of G 77 (m) 
and of the Monte Carlo are entirely correlated and that all other sources of uncertainty are 
uncorrelated. In the energy range where both types of data are available, the ratio of the 
7r°7r + 7r~ to the 77(770) results is constant and equal to 1.14 ± 0.07 (stat.) ± 0.13 (syst.). With 
the improved knowledge of the absolute normalisation provided by the combination of the two 
data sets, the error on the total 77 rate is 10.9%. 
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8 Results 



In this section, the experimental meson rates are presented and compared to the predictions of 
the JETSET and HERWIG Monte Carlo models. To simplify the comparison, the predictions 
of the default versions of JETSET 7.4 and HERWIG 5.9 are used here, except for the a^ which 
is not produced in the default version of JETSET 7.4 and for which the prediction of the 



JETSET version of ref. [19| is used. Two aspects of the model predictions are investigated: the 



shape of the momentum spectra and the integrated rates. 



8.1 Differential cross-sections 



The differential cross-sections as a function of xe and \og(l/x p ) of the photon and the ir°, rj, 
p ± , oj, rj' and ag mesons are obtained by dividing the yields by the corresponding bin widths. 
For the p ± and a^ resonances, the relationship between the meson energy and its momentum 
varies from event to event. The yields are thus evaluated separately in bins of xe and in 
bins of log(l/x p ). The differential cross-sections are listed in tables || to [11] together with 



the statistical and systematic errors. In fig. [T3|, the data are compared to the JETSET 7.4 



predictions normalised to the measured rate. In the simulation, the difference between the 



slopes of the photon and ir° distributions (fig. |l3|a) is strongly constrained by the fact that 92% 



of the photons come from 7r° decays. The measured shapes are consistent with the photon and 
7T° predictions. The p ± and uj mesons, both vector particles with similar masses, have similar 
slopes (fig. |T3"|b) that are also well reproduced by the Monte Carlo. However the predicted r\ 
spectrum is too soft compared to the measurement (fig. |T3|c), while the rf prediction is consistent 
with the data within errors. The slope of the ag distribution (fig. |l3|d) is well reproduced by 
JETSET with the parameters of ref. Jl9| . 



To emphasise the low-momentum portions of the spectra which represent the largest fraction 



of the inclusive rates, the differential cross-sections are presented in fig. as a function of 
log(l/x p ). In this figure, the full and dashed curves are the absolute predictions of default 
JETSET 7.4 and HERWIG 5.9, respectively, except for the a^ prediction of JETSET, taken 
from ref. |1{J. JETSET 7.4 reproduces the photon and 7r° data slightly better than HERWIG 
5.9. The predictions for the a^ are quite similar, but HERWIG 5.9 provides a better description 
of the p ± , uj and rj' data. However, the r) momentum spectra predicted by both models are too 
soft. 



8.2 Maxima of the \og(l/x p ) distributions 

The location of the maximum of the log(l/:Tp) distribution is expected to be correlated with 
the mass of the particle |3"6 |. In addition, its value provides a quantitative measurement of 
the hardness of the momentum spectrum. The value is extracted by fitting a Gaussian to the 
data close to the maximum of the distribution. These values are listed in table O, and shown 



in fig. |15| a where they are compared to other measurements at LEP and to the predictions of 
JETSET and HERWIG. The errors quoted in the table are the sum of the fit errors and of the 
uncertainties due to variations of the range of the Gaussian fits. 
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As can be seen in fig. 15 a,, both JETSET and HERWIG reproduce the data within errors, 
except for the rj, for which the maximum is at a higher value in both JETSET and HERWIG. 
This confirms that the 77 spectrum in the two models is too soft as was already noted in 
section O The OPAL results agree with those of DELPHI @ and L3 PHT21. 



8.3 Rates extrapolated to < xe < 1 



The extrapolation of the rates to the full range of xe is done using the fractions of the rate lying 
outside the measured range predicted by JETSET 7.4, HERWIG 5.9 and a Gaussian fit to the 
log(l/x p ) distributions. The applied corrections correspond to the average of the lowest and 
highest values and the maximum deviation is taken as the systematic error on the procedure. 
The data used to evaluate the extrapolation factors and the associated errors are listed in 



table [13|. The results of the Gaussian fits are considered because of their good description of 
the data within errors and because in at least one case (the 77) it appears that the JETSET and 
HERWIG shapes may not be appropriate. Other shapes that describe the data equally well 
are also considered. For the 7r° and the r\ they are, respectively, the 7r and distributions 
measured by OPAL The extrapolation factors derived from these experimental shapes 

are all within the range covered by the predictions of JETSET, HERWIG and the Gaussian fit. 

The particle multiplicities per hadronic Z° decay extrapolated to the full energy range are: 



(n 7 ) 
(n n o) 

K> 
M 

(ry) 



20.97 ± 0.02 ± 1.07 ± 0.42 
9.55 ± 0.06 ± 0.72 ± 0.21 
0.97 ± 0.03 ± 0.10 ± 0.04 
2.40 ± 0.06 ± 0.43 ± 0.02 
1.04 ± 0.04 ± 0.13 ± 0.03 
0.14 ± 0.01 ± 0.02 ± 0.01 
0.27 ± 0.04 ± 0.10 ± 0.01 



where the first errors are statistical, the second systematic and the third are from the extrap- 
olation procedure. 



8.4 Discussion of the rates and comparison with models 

In table 111], the measured rates are compared to those from other LEP experiments and to the 
predictions of JETSET 7.4 and HERWIG 5.9. The ratio of the measured rates to the JETSET 



7.4 predictions are shown in fig. [T5] b together with the results from other LEP experiments |], 
£J|-T2|]. In this figure, the results obtained by ALEPH in limited energy ranges for the photon, 
7] and rf are divided by the JETSET 7.4 predictions in the corresponding range. All the 
measurements are consistent with each other, except perhaps for the rj' where the measured 
rate agrees with the ALEPH result H, but is more than two standard deviations away from 



that of L3 [12 



According to JETSET 7.4, 97.0% of all photons come from n°, 77, uj and 77' decays. The pre- 
diction of HERWIG 5.9 is 96.0%. In comparison, the sum of the measured 7r°, rj, uj and 7/ rates 
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multiplied by the known photon multiplicities in their decays [^6] accounts for (95 ± 5)% of the 
measured number of photons^. The error on the ratio is calculated assuming that the errors on 
the photon, n° and r\ rates are fully correlated except for those associated with the background 
subtraction and the fit to the invariant mass spectraQ. This good agreement provides further 
evidence that the models offer a reasonable description of the inclusive production of photons 
in Z° decays, an assumption on which the measurement of the photon rates relies. 

The production rates of the ir° and p can be compared with those of their isospin parters. 
Using the measurements for the 7r =1= [fj| and the p° [[],|8j , the following ratios of rates are obtained: 

2(n 7r o)/(n 7r ±) = 1.12 ± 0.01 ± 0.08 , 
2(n p o)/(n p ±) = 1.08 ± 0.04 ± 0.20 , 

where the first errors are statistical and the second systematic. These ratios are consistent with 
the predictions of JETSET 7.4 (1.132 and 1.064, respectively) and HERWIG 5.9 (1.147 and 
1.032). In these models, most of the deviation from unity comes from the decays rj — > 7r°7r 7r° 
and rj' — > p°7. 



As seen in fig. |15|b, the measured rates are consistent with the predictions of JETSET 7.4 
and HERWIG 5.9, except for the u and the rj' for which the rates are more than two standard 
deviations smaller than the JETSET 7.4 prediction. The failure of JETSET to reproduce the 
i]' rate is well understood, as it assumes a similar strange quark content for the rj and rf and 
neglects the effect of their difference in mass on their relative production rate. For this reason, 
the suppression of the rf relative to the rj is a free parameter in JETSET 7.4. The present data 
suggest that the current suppression factor of 0.4 should be further reduced. In contrast to the 
rj', no single parameter can modify the u rate in JETSET independently of all other mesons. 
In that model an increase of the u rate is necessarily accompanied by an equivalent increase of 
the p°, p + and p~ rates because, with ideal mixing, these mesons are the corresponding isospin 
1 = and 1=1 states. Indeed, the measured p and u rates (table [nj) are consistent with 
1 = 0/1=1 symmetry, albeit with a large error. However that symmetry can be broken by 
cascade decays of heavier mesons such as the L = 1 states. The experimental information on 
the production of these states is limited and the present ag data is interesting in this respect. 

Table [14| shows that the a^~ and ag are produced at rates comparable to the rf, a meson 
of equal spin and similar mass. However, the strangeness content of the ag and rf are not 
expected to be the same. A more relevant comparison is with the fo(980) meson, which is the 
isospin 1=0 partner of the ag according to the quark model of mesons. The ratio of the rates 
of the ag and the fg(980) || is 1.9 ± 0.8, compatible with the expected value of 2. 



In the HERWIG cluster fragmentation model pi] , the relative production rates of light- 
flavour mesons are mostly determined by their masses, which affects the phase space available 
for the cluster decay. This simple ansatz appears to be able to reproduce the measured ag 
rate. The a^ is not present in the default version of JETSET. The inclusion of L = 1 meson 
production in hadronisation requires the tuning of additional parameters. The predictions 
shown in fig. [15] correspond to the choice of parameters of ref. |19| optimised, in part, to 



6 In the data, the 7r° and rj decays alone account for approximately 91% and 4% of all observed photons, 
respectively, in good agreement with the JETSET 7.4 predictions of 91.9% and 4.0%. 

7 These are the errors labelled as "Background subtraction" , "Gaussian peaks" and "Background normalisa- 
tion range" in tables 1 to 3. The contributions of the ui and the rf to the total error are negligible. 
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reproduce the available data on D and B mesons. The agreement for the ag may be accidental 
since the parameters of ref. Jl9| also predict a substantial bi(1232) ±,C) rate of 0.92 meson per 



Z° decays. The present u and p data do not support this prediction. With the bi(1232) ±,ci 
decaying exclusively to ujtt, the parameters of ref. [0 result in an u rate which exceeds that of 
its 7 = 1 partner, the p, with (nj) — (n p ±)/2 = 0.56. The measured difference is —0.17 ± 0.26, 
in better agreement with the prediction of —0.06 of default JETSET 7.4 with no L = 1 mesons, 
and with the prediction of 0.00 of HERWIG 5.9, with a bi(1232) ± '° rate of 0.32. Thus, despite 
the lack of direct measurements for several L = 1 meson states, it appears possible to constrain 
the JETSET model by using the available data on all other mesons. 



9 Conclusion 



The inclusive particle multiplicity per hadronic Z° decay and the differential cross-section have 
been measured for photons and for 7r°, 77, p ± , to, rf and a^ mesons. The a^ is observed for 
the first time in high-energy e + e~ collisions. It is produced at a rate comparable to that of 
mesons with a similar mass, such as the rf and the /o(980). The inclusive p production is 
measured for the first time in hadronic Z° decays. The models JETSET 7.4 and HERWIG 5.9 
with their default parameters reproduce the shape of the measured differential cross-sections, 
with the exception of that of the r\ meson which is too soft in both models. The absolute rates 
in HERWIG 5.9 are in good agreement. In JETSET 7.4, the production rates of the u and 77 
are overestimated by 20% and 50%, respectively. The present ag data is a valuable input for 
the determination of the parameters required for the inclusion of L = 1 mesons in JETSET. 
These parameters are further constrained by the data on other mesons like the u and the p ± . 
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Photon sample 


7 


7c 


7 + 7c 


xe range 


0.003 - 1.000 


0.003 - 1.000 


0.003 - 1.000 


Integrated rate 


16.79 


16.96 


16.84 


Errors (%) 




Statistics (data) 


0.1 


0.2 


0.1 


Statistics (Monte Carlo) 


0.1 


0.1 


0.1 


Difference between Monte Carlos 


3.4 


3.0 


3.3 


P-y(S) variations 


0.8 




0.6 


Background subtraction 


0.8 


0.5 


0.7 


Nuclear interactions 


2.0 


0.9 


1.7 


Energy scale 


1.5 


1.3 


1.4 


Track and cal. simulation 


1.8 


7.6 


3.1 


Total error (%) 


4.7 


8.3 


5.1 



Table 1: Number of photons per hadronic Z decay in the xe range covered by the measurement 
together with its statistical and systematic uncertainties (in %). The three columns are the 
results obtained with the calorimetric sample (7), the conversion sample (7J and the combined 
sample (7 + 7 C ). 
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7T° sample 


77 


77c 


7c7c 


77 + 77c + 7c7c 


xe Range 


0.007 - 0.500 


0.007 - 0.500 


0.009 - 0.300 


0.007-0.500 


Integrated rate 


8.37 


8.12 


7.51 


8.29 


Rate (0.009 < x E < 0.3) 


7.80 


7.35 


7.51 


7.65 


Errors (%) 




Statistics (data) 




0.6 


2.1 


4.9 


0.6 


Statistics (Monte Carlo) 




0.4 


1.0 


2.8 


0.4 


Difference between Monte 


Carlos 


4.7 


4.0 


5.9 


4.1 


P 77 variations 




2.0 


2.7 


9.3 


2.2 


Background subtraction 




3.7 


7.2 


11.5 


3.8 


Gaussian peaks 




0.8 


1.1 


3.0 


0.8 


Nuclear interactions 




1.6 


0.7 


0.1 


1.3 


Energy scale (1%) 




1.8 


1.9 


1.5 


1.6 


Background normalisation 


range 


4.3 


2.3 


5.0 


3.1 


Track and cal. simulation 




2.2 


7.6 


15.2 


2.5 


Total error (%) 


8.4 


12.2 


23.6 


7.6 



Table 2: Statistical and systematic uncertainties (in %) on the number of 7r° mesons per 
hadronic Z° decay measured using the individual 7r° — > 77, n° — > 77 c and 7r° — > 7 C 7 C samples, 
and the combined sample. 



i] sample 


77 


77c 


77 + 77c 


xe Range 


0.025 - 1.000 


0.040 - 0.300 


0.025 - 1.000 


Integrated rate 


0.746 


0.569 


0.745 


Rate (0.04 <x E < 0.30) 


0.558 


0.569 


0.558 


Errors (%) 




Statistics (data) 


2.6 


6.5 


2.6 


Statistics (Monte Carlo) 


1.7 


4.5 


1.7 


Difference between Monte Carlos 


3.9 


9.7 


3.8 


P 77 variations 


10.5 


26.0 


10.5 


Background subtraction 


7.4 


16.7 


7.3 


Gaussian peaks 


0.2 


1.1 


0.2 


Energy scale (1%) 


0.8 


0.4 


0.8 


Background normalisation range 


1.9 


3.9 


1.9 


Track and cal. simulation 


1.1 


6.9 


1.1 


Total Error (%) 


14.0 


34.3 


13.9 



Table 3: Statistical and systematic uncertainties (in %) on the number of 77 mesons per 
hadronic Z° decay measured using the individual rj — > 77 and r\ — > 77,, samples, and the 
combined sample. 
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Meson 


V 


P ± 




rf 


aj 


Decay mode 


7r°7r + 7r~ 


ttV 








xe range 


0.025-0.4 




0.025-0.6 


0.05-0.8 




log(l/xp) range 




0.0-5.0 






0.0-3.5 


Rate 


0.898 


2.36 


0.883 


0.103 


0.214 


Errors (%) 




Statistics (data) 


6.5 


2.3 


3.6 


10.1 


16.8 


Statistics (Monte Carlo) 


1.7 


1.7 


0.5 


1.6 


2.5 


Difference between Monte Carlos 


3.7 




5.9 


4.9 




G 77 (m) variations 


2.4 


10.9 


3.9 


7.2 


13.6 


Background subtraction 


4.5 


7.4 


3.0 


5.0 


30.4 


Mass resolution 


4.0 


7.0 


8.1 


6.8 




Breit-Wigner extrapolation 




3.1 








Track and cal. simulation 


4.0 


4.0 


4.0 


3.0 


3.0 


Partial reconstruction 


6.6 


9.0 


5.8 




6.6 


Branching Ratio 


2.5 






3.4 


10.3 


Total Error (%) 


12.9 


18.4 


13.7 


L 16.5 


L 39.4 



Table 4: Statistical and systematic uncertainties (in %) on the number of r), p ± , to, rf and a 
mesons per hadronic Z° decay measured with the channels combining a ir° or an 77 meson with 
charged pions. For the p ± and ag resonances, the error associated with the difference between 
the Monte Carlos is included in the error associated with the background subtraction. 
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x E 


I da 


ln(l/xp) 


1 da 


Vhad dXft 


a had d\n(l/xp) 


o nns - 

u.uuo 


U.UUt: 


1309 ± 13 ± 330 


O .O X 


O.O^ 


4.55 ± 0.05 ± 1.15 


n nn4 - 

U.UUt: 


007 

U.UU I 


986 ± 2 ± 78 


o .o^ 


4 Qfi 

t:. C/U 


5.29 ± 0.01 ± 0.42 


n no7 

U.UU I 


00Q 


749 ± 1 ± 34 


4 Qfi 


4 71 


5.96 ± 0.01 ± 0.27 


00Q 

U.UUc/ 


011 


613 ± 1 ± 28 


4 71 


4 ^1 


6.11 ± 0.01 ± 0.28 


01 1 

U .U ± ± 


01 ^ 

U.UIO 


508 ± 1 ± 23 


4 ^1 


4 ^4 


6.08 ± 0.02 ± 0.27 


01^ 


01 

U.UIU 


404 ± 1 ± 16 


4 ^4 


4 14 


5.84 ± 0.01 ± 0.23 


01 (\ 

U .U 1U 


090 
u .uz,u 


303 ± 1 ± 11 


4 14 


3 Q1 


5.43 ± 0.01 ± 0.20 


090 
u.uz,u - 


09^ 


225 ± 1 ± 9 


3 Q1 

O . a 1 




5.03 ± 0.01 ± 0.19 


09^ 


0^0 
u.uou 


167 ± 1 ± 7 


o.uy 


3 ^1 

O. O -L 


4.58 ± 0.01 ± 0.18 


0^0 


0^ 
u.uoo 


131 ± 1 ± 5 


^ ^1 
o.oi - 


o.oo 


4.24 ± 0.01 ± 0.15 


0^ 


040 

U.Ut:U 


103 ± 1 ± 4 


o.oo 


3 99 


3.87 ± 0.01 ± 0.15 


040 - 

U.UtiU 


0^0 


76.9 ± 0.2 ± 2.8 


^ 99 - 


3 00 

O.UU 


3.45 ± 0.01 ± 0.13 


0^0 
u.uou 


OfiO 
u.uuu 


53.4 ± 0.2 ± 2.0 


3 00 
o.uu 


9 81 


2.93 ± 0.01 ± 0.11 


OfiO 


070 
U.U i u 


39.2 ± 0.2 ± 1.4 


9 81 




2.54 ± 0.01 ± 0.09 


070 


08^ 
u.uoo 


28.1 ± 0.1 ± 1.0 




9 47 


2.17 ± 0.01 ± 0.08 


08^ 


100 

U . 1UU 


19.6 ± 0.1 ± 0.8 


9 47 


9 ^0 


1.81 ± 0.01 ± 0.08 


1 00 


1 9^ 


13.0 ± 0.1 ± 0.5 


9 ^0 


9 08 


1.45 ± 0.01 ± 0.06 


1 9^ 


1 ^0 

U . J-OU 


8.04 ± 0.05 ± 0.35 


9 OR 

Zi.UO 


1 QO 

1 . c/U 


1.103 ± 0.006 ± 0.048 


0.150 - 


0.200 


4.50 ± 0.03 ± 0.26 


1.90 - 


1.61 


0.783 ± 0.005 ± 0.045 


0.200 - 


0.300 


1.71 ± 0.02 ± 0.19 


1.61 - 


1.20 


0.422 ± 0.004 ± 0.047 


0.300 - 


0.400 


0.507 ± 0.010 ± 0.087 


1.20 - 


0.92 


0.176 ± 0.003 ± 0.030 


0.400 - 


0.500 


0.184 ± 0.005 ± 0.038 


0.92 - 


0.69 


0.082 ± 0.002 ± 0.017 


0.500 - 


0.600 


0.065 ± 0.002 ± 0.011 


0.69 - 


0.51 


0.036 ± 0.001 ± 0.006 


0.600 - 


0.800 


0.017 ± 0.001 ± 0.002 


0.51 - 


0.22 


0.012 ± 0.000 ± 0.002 


0.800 - 


1.000 


0.0023 ± 0.0003 ± 0.0010 


0.22 - 


0.00 


0.0020 ± 0.0003 ± 0.0009 



Table 5: Photon fragmentation function obtained by combining the calorimeter and conversion 
results. The quoted errors are statistical and systematic, respectively. 
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x E 


1 da 


ln(l/x p ) 


1 




dc 


r 


a had dXE 


ahad <fln(l/x p ) 


0.007 - 


0.009 


254 + 1 8 + 48 


5.06 - 


A 1-71-7 

4.77 


1 >T A I 

1.74 ± 





12 


± 0.33 


0.009 - 


0.011 


266 ± 12 ± 38 


4.77 - 


4.55 


2.42 ± 





11 


± 0.34 


0.011 - 


0.013 


248 + 6 + 28 

O _l \j _l aiU 


4.55 - 


4.37 


2.78 ± 





06 


± 0.31 


0.013 - 


0.016 


21 1 +3 + 18 

JL 1. _l '.J _l X O 


4.37 - 


4.15 


2.92 ± 





05 


± 0.25 


0.016 - 


0.020 


i 70 _i_ o + 1 4 


4.15 - 


3.92 


3. II ± 





03 


± 0.25 


0.020 - 


0.025 


139 ±2 + 6 


3.92 - 


3.70 


3.05 ± 





03 


± 0.13 


0.025 - 


0.030 


113 ±1 + 5 


3.70 - 


3.51 


3.06 ± 





03 


± 0.13 


0.030 - 


0.035 


Q4 1 ± 9 ± 4 


3.51 - 


3.36 


3.03 ± 





03 


± 0.13 


0.035 - 


0.040 


77 7 + 08 + 43 


3.36 - 


3.22 


2.89 ± 





03 


± 0.16 


0.040 - 


0.050 


62 5 + 4 + 3 9 


3.22 - 


3.00 


2.79 ± 





02 


± 0.17 


0.050 - 


0.060 


45 7 ± 4 ± 3 


3.00 - 


2.81 


2.50 ± 





02 


± 0.16 


0.060 - 


0.070 


34 7 + 3 + 3 


2.81 - 


2.66 


2.25 ± 





02 


± 0.20 


070 - 

w . w i yj 


08^ 


26 2 + 2 + 1 8 


2 66 - 

Zj . \J\J 


2 47 


2.02 ± 





02 


± 0.14 


0.085 - 


0.100 


19.4 ± 0.2 ± 1.4 


2.47 - 


2.30 


1.79 ± 





02 


± 0.13 


0.100 - 


0.125 


13.2 ± 0.1 ± 2.9 


2.30 - 


2.08 


1.48 ± 





02 


± 0.32 


0.125 - 


0.150 


9.05 ± 0.13 ± 0.76 


2.08 - 


1.90 


1.240 ± 0.017 ± 0.105 


0.150 - 


0.200 


5.36 ± 0.10 ± 0.69 


1.90 - 


1.61 


0.931 ± 0.017 ± 0.120 


0.200 - 


0.300 


2.26 ± 0.13 ± 0.38 


1.61 - 


1.20 


0.558 ± 0.031 ± 0.094 


0.300 - 


0.400 


0.764 ± 0.085 ± 0.309 


1.20 - 


0.92 


0.266 ± 0.030 ± 0.107 


0.400 - 


0.500 


0.455 ± 0.095 ± 0.244 


0.92 - 


0.69 


0.204 ± 0.043 ± 0.110 



Table 6: ir° fragmentation function obtained by combining the 7r° — > 77, 7r° — > 77 c and 
7T° — > 7 C 7 C data. The quoted errors are statistical and systematic, respectively. 



x E 


1 da 


ln(l/ 


Xp) 




1 


da 






ahad dxp 




ahad cfln(l/ 






0.025 - 


0.035 


10.6 ± 1.5 ± 2.4 


3.82 - 


3.42 


0.261 


± 


0.038 


± 


0.046 


0.035 - 


0.050 


7.63 ± 0.78 ± 1.27 


3.42 - 


3.03 


0.294 


± 


0.030 


± 


0.038 


0.050 - 


0.075 


5.10 ± 0.38 ± 0.61 


3.03 - 


2.60 


0.302 


± 


0.023 


± 


0.028 


0.075 - 


0.100 


3.81 ± 0.21 ± 0.44 


2.60 - 


2.31 


0.324 


± 


0.018 


± 


0.032 


0.100 - 


0.125 


2.83 ± 0.12 ± 0.28 


2.31 - 


2.08 


0.314 


± 


0.014 


± 


0.028 


0.125 - 


0.150 


2.21 ± 0.10 ± 0.22 


2.08 - 


1.90 


0.301 


± 


0.014 


± 


0.027 


0.150 - 


0.200 


1.46 ± 0.05 ± 0.13 


1.90 - 


1.61 


0.252 


± 


0.009 


± 


0.021 


0.200 - 


0.300 


0.733 ± 0.026 ± 0.062 


1.61 - 


1.20 


0.180 


± 


0.006 


± 


0.014 


0.300 - 


0.400 


0.364 ± 0.022 ± 0.047 


1.20 - 


0.92 


0.126 


± 


0.008 


± 


0.014 


0.400 - 


0.500 


0.220 ± 0.019 ± 0.031 


0.92 - 


0.69 


0.099 


± 


0.008 


± 


0.011 


0.500 - 


0.600 


0.086 ± 0.010 ± 0.019 


0.69 - 


0.51 


0.047 


± 


0.006 


± 


0.009 


0.600 - 


0.800 


0.033 ± 0.004 ± 0.008 


0.51 - 


0.22 


0.023 


± 


0.003 


± 


0.005 


0.800 - 


1.000 


0.0013 ± 0.0004 ± 0.0011 


0.22 - 


0.00 


0.0012 ± 0.0004 ± 0.0009 



Table 7: r] fragmentation function obtained by combining the 77, 77 c and 7r°7r + 7r data. The 
quoted errors are statistical and systematic, respectively. 
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X E 




1 da 


ln(l/ 






l 


da 










° had dx p 




"had cfln(l/ 






n 


016 - 





025 


17 3 + 81 +122 


5 0- 

xj . W 


4 5 


1 71 


± 


0.008 




081 


n 


025 - 





035 


32 3 + 25 + Q7 


4 5- 


4 


41 9 


± 


0.035 


4 


111 


n 


035 - 





050 


oi q + n 7 + 4 5 


4 0- 


3 5 


500 


± 


0.092 


4 


1 38 


n 


D5D - 

\JO\J 


n 

u 


075 


i a 7 _i_ n a _i_ i c 


3 5- 


3 


692 


4- 

1 


fl 028 


-|- 


165 





075 - 





100 


9.89 ± 0.40 ± 1.46 


3.0 - 


2.5 


0.868 


± 


0.021 


± 


0.126 





100 - 





125 


7.11 ± 0.25 ± 1.04 


2.5 - 


2.0 


0.805 


± 


0.022 


± 


0.104 





125 - 





150 


5.90 ± 0.25 ± 0.78 


2.0 - 


1.5 


0.603 


± 


0.017 


± 


0.078 





150 - 





200 


3.60 ± 0.12 ± 0.48 


1.5 - 


1.0 


0.419 


± 


0.014 


± 


0.073 





200 - 





300 


2.02 ± 0.07 ± 0.21 


1.0 - 


0.5 


0.217 


± 


0.010 


± 


0.055 





300 - 





400 


1.03 ± 0.04 ± 0.27 


0.5 - 


0.0 


0.034 


± 


0.004 


± 


0.019 





400 - 





600 


0.430 ± 0.023 ± 0.081 



















600 - 





800 


0.075 ± 0.013 ± 0.032 



















800 - 


1 


000 


0.013 ± 0.003 ± 0.009 

















Table 8: p fragmentation function. The quoted errors are statistical and systematic, re- 
spectively. Because of the width of the p ± , the relation between xe and x p varies with mass. 
Therefore the analysis is performed first with bins of xe (first two columns) and then repeated 
with bins of \og{l/x p ). 



x E 


1 da 


ln(l/x p ) 




i 


da 






&had dx p 




<?had <21n(l/ 


<->) 




0.025 - 


0.035 


15.2 ± 2.4 ± 2.1 


4.01 - 


3.49 


0.293 


± 


0.046 


± 


0.040 


0.035 - 


0.050 


9.88 ± 0.84 ± 1.48 


3.49 - 


3.06 


0.344 


± 


0.029 


± 


0.051 


0.050 - 


0.075 


5.82 ± 0.35 ± 0.75 


3.06 - 


2.62 


0.330 


± 


0.020 


± 


0.043 


0.075 - 


0.100 


4.12 ± 0.25 ± 0.54 


2.62 - 


2.32 


0.344 


± 


0.021 


± 


0.045 


0.100 - 


0.125 


2.74 ± 0.16 ± 0.32 


2.32 - 


2.09 


0.299 


± 


0.018 


± 


0.035 


0.125 - 


0.150 


2.23 ± 0.14 ± 0.24 


2.09 - 


1.90 


0.301 


± 


0.018 


± 


0.032 


0.150 - 


0.200 


1.45 ± 0.09 ± 0.17 


1.90 - 


1.61 


0.250 


± 


0.016 


± 


0.029 


0.200 - 


0.300 


0.789 ± 0.049 ± 0.099 


1.61 - 


1.21 


0.193 


± 


0.012 


± 


0.024 


0.300 - 


0.400 


0.335 ± 0.037 ± 0.042 


1.21 - 


0.92 


0.116 


± 


0.013 


± 


0.014 


0.400 - 


0.600 


0.130 ± 0.027 ± 0.028 


0.92 - 


0.51 


0.064 


± 


0.013 


± 


0.014 



Table 9: uo fragmentation function. The quoted errors are statistical and systematic, respec- 
tively. 
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x E 



0.050 
0.070 
0.100 
0.125 
0.150 
0.200 
0.300 
0.400 
0.600 



0.070 
0.100 
0.125 
0.150 
0.200 
0.300 
0.400 
0.600 
0.800 



1 da 

Zhad dx p 



1.01 ± 
0.462 ± 
0.460 ± 
0.293 ± 
0.354 ± 
0.137 ± 
0.088 ± 
0.034 ± 
0.013 ± 



0.38 ± 0.14 
0.180 ± 0.073 
0.144 ± 0.082 
0.099 ± 0.049 
0.068 ± 0.054 
0.028 ± 0.017 
0.020 ± 0.011 
0.010 ± 0.006 
0.006 ± 0.003 



\n(l/x p ) 



3.09 
2.71 
2.33 
2.09 
1.91 
1.61 
1.21 
0.92 
0.51 



2.71 
2.33 
2.09 
1.91 
1.61 
1.21 
0.92 
0.51 
0.22 



da 



^had dln(l/x p ) 



0.052 
0.036 
0.050 
0.039 
0.061 
0.034 
0.030 
0.017 
0.009 



± 0.020 ± 
± 0.014 ± 
± 0.016 ± 
± 0.013 ± 
± 0.012 ± 
± 0.007 ± 
± 0.007 ± 
± 0.005 ± 
± 0.004 ± 



0.007 
0.006 
0.009 
0.007 
0.009 
0.004 
0.004 
0.003 
0.002 



Table 10: rj' fragmentation function. The quoted errors are statistical and systematic, respec- 
tively. 



X E 


1 da 


ln(l/x p ) 


1 da 


Chad dx p 


ahad d\n(l/x p ) 


0.050 - 0.070 
0.070 - 0.100 
0.100 - 0.125 
0.125 - 0.150 
0.150 - 0.200 
0.200 - 0.300 
0.300 - 0.400 
0.400 - 0.600 
0.600 - 0.800 
0.800 - 1.000 


1.65 ± 1.03 ± 0.75 
1.05 ± 0.49 ± 0.73 
0.747 ± 0.215 ± 0.214 
0.985 ± 0.238 ± 0.560 
0.623 ± 0.107 ± 0.171 
0.207 ± 0.046 ± 0.069 
0.093 ± 0.027 ± 0.040 
0.038 ± 0.015 ± 0.015 
0.014 ± 0.005 ± 0.006 
0.0040 ± 0.0018 ± 0.0024 


3.50 - 3.00 
3.00 - 2.50 
2.50 - 2.00 
2.00 - 1.50 
1.50 - 1.00 
1.00 - 0.50 
0.50 - 0.00 


0.093 ± 0.063 ± 0.050 
0.104 ± 0.041 ± 0.041 
0.076 ± 0.019 ± 0.030 
0.088 ± 0.013 ± 0.023 
0.040 ± 0.009 ± 0.012 
0.019 ± 0.006 ± 0.007 
0.0071 ± 0.0025 ± 0.0022 



Table 11: ag fragmentation function. The quoted errors are statistical and systematic, re- 
spectively. Because of the width of the aQ , the relation between xe and x v varies with mass. 
Therefore the analysis is performed first with bins of xe (first two columns) and then repeated 
with bins of log(l/:Tp). 
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Particle 


Location of the maximum of da/dln(l/x p ) 


OPAL 


L3 [g^g 


DELPHI 


El 


JETSET 7.4 


HERWIG 5.9 


7 
V 

p ± 

OJ 
7]' 

a 


4.61 ± 0.12 
3.77 ± 0.11 
2.64 ± 0.14 
2.63 ± 0.15 
2.89 ± 0.24 
2.21 ± 0.42 
2.57 ± 0.50 


3.90 ±8;fl 
2.52 ± 0.12 

2.86 ± 0.20* 
2.69 ± 0.10* 


3.96 ± 0.13 


4.54 
3.78 
2.94 
2.69 
2.77 
2.48 
2.62 


4.63 
3.86 
3.01 
2.70 
2.77 
2.17 
2.72 



Table 12: Location of the maximum of the log(l/x p ) distributions, determined from a Gaussian 
fit to the data in the region around the maximum. The OPAL measurements are compared 
to other LEP measurements and to the predictions of the default versions of JETSET 7.4 and 
HERWIG 5.9. The errors on the predictions are typically ±0.01. The results marked with 
an asterisk are extracted assuming that the shape of the log(l/x p ) distribution is given by a 



MLLA calculation 12 





7 


7T° 


V 




10 


ri 


a o = 


Measured range: 
















Min. log(l/a;p) 


0.00 


0.69 


0.00 


0.00 


0.51 


0.22 


0.00 


Max. \og{l/x p ) 


5.81 


5.06 


3.82 


5.00 


4.01 


3.09 


3.50 


% of rate in range 
















JETSET 7.4 


81.9 


88.7 


78.7 


98.9 


87.7 


68.9 


80.5* 


HERWIG 5.9 


79.4 


87.7 


77.9 


98.7 


87.1 


74.5 


78.6 


Gaussian fit 


78.7 


84.9 


82.2 


97.4 


82.3 


77.7 


81.6 


7r =1= shape 




88.4 












Kg shape 






82.2 










Combined % 


80.3 


86.8 


81.1 


98.2 


85.0 


73.3 


80.1 


Error 


1.6 


1.9 


3.1 


0.8 


2.7 


4.4 


1.5 


Measured rate 


16.84 


8.29 


0.789 


2.36 


0.883 


0.103 


0.214 


Extrapolated rate 


20.97 


9.55 


0.973 


2.40 


1.039 


0.141 


0.267 


Extrapolation error 


0.42 


0.21 


0.038 


0.02 


0.033 


0.008 


0.005 



Table 13: Data used for the extrapolation to the unobserved energy/momentum ranges. The 
result marked with an asterisk is from the JETSET tune of ref . |L9| . 
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Experimental results 


JETSET 
7.4 


HERWIG 
5.9 


OPAL 


ALEPH § 


DELPHI [|] 


L3 


photon 
xe range 
N~ in range 
N 7 all xe 


0.003-1.000 
16.84 ± 0.86 
20.97 ± 1.15 


0.018-0.450 
7.37 ± 0.24 






20.76 


22.65 


7T° 

xe range 
N n o in range 
N w o all ie 


0.007-0.400 
8.29 ± 0.63 
9.55 ± 0.76 


0.025-1.000 
4.80 ± 0.32 
9.63 ± 0.64 


0.011-0.750 

7.1 ± 0.8 

9.2 ± 1.0 


0.004-0.150 
8.38 ± 0.67 
9.18 ± 0.73 


9.60 


10.29 


V 

xe range 
AL in ranee 

all xe 
N v x p > 0.1 


0.025-1.000 
0.79 ± 0.08 
0.97 ± 0.11 
0.344 ± 0.030 


0.100-1.000 
0.282 ± 0.022 

0.282 ± 0.022 




0.020-0.300 
0.70 ± 0.08 
0.91 ± 0.11 


1.00 
0.286 


0.92 
0.243 


P± 
xe range 
AT ± in range 
A/p± all xe 


0.016-1.000 
2.36 ± 0.42 
2.40 ± 0.44 








2.82 


2.29 


uj 

xe range 
AT W in range 
N w all 


0.025-0.800 
0.88 ± 0.12 
1.04 ± 0.14 


0.053-1.000 
0.64 ± 0.08 
1.07 ± 0.14 




0.026-0.300 
0.94 ± 0.14 
1.17 ± 0.17 


1.35 


1.14 


rf 

xe range 
A^/ in range 

N v i all 
iV; x p > 0.1 


0.050-0.800 
0.103 ± 0.017 

0.14 ± 0.03 
0.069 ± 0.012 


0.100-1.000 
0.064 ± 0.014 

0.064 ± 0.014 




0.023-0.240 
0.25 ± 0.04 


0.297 
0.127 


0.122 
0.060 


a o 

i£ range 
N+ in range 
N+ all xe 


0.050-1.000 
0.21 ± 0.08 
0.27 ± 0.11 








0.210 


0.221 



Table 14: Summary of the measurements of the particles rates described in this paper compared 
to other measurements at LEP and to JETSET and HERWIG predictions. 
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Figure 1: Distribution of variables used for the discrimination of photons recorded in the 
electromagnetic calorimeter. The points represent the data and the histogram represents the 
Monte Carlo simulation, normalised to one event. The Monte Carlo photons and background 
are shown as hatched and dashed histograms, respectively. The insert shows the weights used 
in eq. [3]. The variables are: a) the photon energy; b) the photon energy multiplied by the angle 
to the closest charged track; c) the minimum of the invariant mass of the photon with any other 
photon; d) the fitted shower energy divided by the sum of the energy in the 3x3 array of blocks 
around the shower; e) the distance of the shower relative to the centre of the block, a value of 
1 corresponding to the edge of the block; f) the maximum value of (E-i — E )(E + i — E )/Eq, 
where E_ ly E , E +1 are the energies deposited in 3 consecutive blocks in either 9 or 0, the 
index corresponding to the block where the shower is centered; and g) the sum of the residual 
of the shower fit in the 3x3 array of blocks around the shower divided by the shower energy. 
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Figure 2: Distribution of P 7 , the variable used for the selection of calorimeter photons in 
the measurement of the photon rates. The points are data and the histogram represents the 
Monte Carlo simulation. The Monte Carlo signal (dashed histogram) and background (hatched 
histogram) are also shown. The Monte Carlo signal has been normalised to the measured photon 
rate and the background to the number of events in the data sample. 
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Figure 3: Radial coordinate r of the photon conversion, in the data (dots) and in the Monte 
Carlo (histogram). The shaded histogram shows the contribution from background, according 
to the Monte Carlo. The Monte Carlo sample is normalised to the same number of events. The 
peaks close to 9, 18 and 23 cm correspond to concentration of material in the detector. The 
abrupt cut at r = 50 cm is due to an additional cut on the dE/dx of the two electron tracks. 
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Figure 4: a-f: Distribution of P 77 for the channels 7r° — > 77, 77 — > 77, 7r° — > 77 c , 77 — > 77 c , 
7T° — > 7 C 7 C and r] — > 7 C 7 C . The points represent the data and the histograms represent the 
Monte Carlo simulation normalised to the same number of events. The Monte Carlo signals 
and backgrounds are shown separately as hatched and dashed histograms, respectively. 
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Figure 5: Absolute efficiency as a function of the scaled meson energy Xe for the reconstruction 
and identification of n° and r\ mesons, according to the Monte Carlo simulation. The efficiencies 
for the decay channels 7r° — > 77, 7T° — > 77 c , 7r° — > 7 C 7 C , 77 — > 77 and — > 77 c are shown 
separately. A cut on P 77 > 0.1 is applied for the 7r°, and P 77 > 0.05 for the 77 (see text). 
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Figure 6: a-f) Fits to the invariant mass distribution of pairs of photons, for the decay 
channels 7r° — > 77, 7r° — > 77 c and 7r° — > 7 C 7 C for the lowest and highest P 77 cut used in the 
present analysis. The points represent the OPAL data and the full lines the fits to the data 
where the background (dashed lines) is parameterized as a second-order polynomial. The dotted 
lines correspond to the background evaluated using the Monte Carlo distributions, which may 
contain a small fraction of the signal in cases where the association between photon candidates 
and true photons is ambiguous (see text). 
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Figure 7: a-d) Fits to the invariant mass distribution of pairs of photons, for the decay channels 
7T° — > 77 and i] — > 77 c for the lowest and highest P 77 cut used in the present analysis. The 
points represent the OPAL data and the full lines the fits to the data where the background 
(dashed lines) is parameterized as a second-order polynomial plus a Gaussian for the reflection 
from uo — > 77T° decays. The dotted lines correspond to the background evaluated using the 
Monte Carlo distributions, which may contain a small fraction of the signal in cases where the 
association between photon candidates and true photons is ambiguous (see text). 
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Figure 8: a-d) Mass peak for the decay channels 7r° — > 77, 7r° — > 77 c , 7r° — > 7 C 7 C and 77 — > 77. 
The distributions are obtained by subtracting the fitted background shape from the raw spectra. 
The points represent the OPAL data and the full curves the fit to the Monte Carlo. The Monte 
Carlo distributions are normalised to the area in the data. A cut on P 77 > 0.3 is applied for the 
7T°, and P 77 > 0.2 for the 77 (see text). 
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Figure 9: Fits to the invariant mass distribution of combinations for two different cuts 

on G 77 (m). The points represent the OPAL data and the full curves the background obtained 
in the fit to the data. The dashed histograms correspond to the background evaluated using the 
Monte Carlo distributions. In the lower part of the figures, the signal obtained by subtracting 
the average of the two background from the data is shown as the open circles, and the signal 
in the Monte Carlo, normalised to the same number of events, is shown as the full histogram. 
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Figure 10: Fits to the invariant mass distribution of r/n^ 1 combinations for two different cuts 
on Gr 77 (m). The points represent the OPAL data and the full curves the background obtained 
in the fit to the data. The dashed histograms correspond to the background evaluated using the 
Monte Carlo distributions. In the lower part of the figures, the signal obtained by subtracting 
the average of the two background from the data is shown as the open circles and is multiplied 
by a factor 3 for clarity. The full histogram is the Monte Carlo signal normalised to the same 
number of events. 
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Figure 11: Fits to the invariant mass distribution of 7r°7r + 7r~ combinations, for a cut of 
G 77 (m)> 0.3. The points represent the OPAL data and the full curves the fit to the data. The 
dashed histograms correspond to the background evaluated using the Monte Carlo distributions. 
Fig. (a) and (c): fits in the region of the uj and the rj, respectively. The signal extracted from 
the fit to the data, multiplied by a factor 5 for clarity, is shown by the open circles. Fig. (b) 
and (d): the points are the components proportional to A (see section |6.4| ) extracted from the 
data in the region of the uj and the 77, respectively. The full curves are the fits to the data. 
The dashed histogram is the simulated background, including partially reconstructed mesons 
which appear as a bump close to the signal peak. The dotted curves are the contributions from 
partial reconstructions included in the fit. 
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Figure 12: Fit to the invariant mass distribution of i]7t + tt~ combinations in the data (left) and 
in the Monte Carlo (right). A cut on G 77 (m)> 0.1 has been applied. The curves represent the 
fit and the Monte Carlo background is shown as a dashed histogram. In the lower part of the 
figures, the open symbols represent the signal extracted from the fit. 
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Figure 13: Differential cross-section as a function of the scaled energy xe for the production in 
hadronic Z° decays of (a) photons and tt° mesons, (b) p^ and u mesons, (c) r] and rf mesons and 
(d) ag mesons. The points are the data and the curves are the predictions of default JETSET 
7.4, with a normalisation fitted to the data. The error bars represent the quadratic sum of the 
statistical and systematic errors. 
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Figure 14: OPAL measurement of the inclusive production of photons, 7r° i], uj, rj' and 
a^ mesons in hadronic Z° decays, as a function of \og(l/x p ). The full and dashed lines are the 
absolute predictions of JETSET 7.4 and HERWIG 5.9, respectively, without any additional 
normalisation. The error bars represent the quadratic sum of the statistical and systematic 
errors. 
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Figure 15: a) Position of the maximum of the log(l/:Ep) distributions. The black points are 
the OPAL measurements, and the full and dashed lines are the predictions of JETSET 7.4 and 
HERWIG 5.9, respectively. The measurements of DELPHI |J and L3 []T0|— 12[| are shown as 



white squares and circles, respectively. The maxima are obtained by a Gaussian fit to the data 
close to the maximum, except for L3, which makes model-dependent assumptions concerning 
the shape of the distribution, b) Particle multiplicities (extrapolated to < xe < 1), divided 
by the prediction of JETSET 7.4. The full circles represent the OPAL measurements, and the 
open circles, triangles and squares those of other LEP experiments P, j9[-|12l. The dotted line 
represents the prediction of HERWIG 5.9. 
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